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We analyze the equilibrium frequency-dependent spin current noise and spin conductance through 
a quantum dot in the local moment regime. Spin current correlations behave markedly differently 
from charge correlations. Equilibrium spin correlations are characterized by two universal scal- 
ing functions in the absence of an external field: one of them is related to charge correlations, 
while the other one describes cross-spin correlations. We characterize these functions using a com- 
bination of perturbative and non-pert urbative methods. We find that at low temperatures spin 
cross-correlations are suppressed at frequencies below the Kondo scale, Tk, and a dynamical spin 
accumulation resonance is found at the Kondo energy, uj ^ Tk- At higher temperatures, T > Tr, 
surprising low-frequency anomalies related to overall spin conservation appear in the spin noise and 
spin conductance, and the Korringa rate is shown to play a distinguished role. The transient spin 
current response also displays universal and singular properties. 

PACS numbers: 72.25.-b, 73.63.Kv, 72.15.Qm, 72.70.+m 
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I. INTRODUCTION 

Nanoelectronic devices are likely to provide our future 
technology and serve as basic tools for storing informa- 
tion, quantum computation^'^ or spin manipulation^!^ 
Due to recent developments in fabrication, it is now pos- 
sible to produce and measure spin accumulation in meso- 
scopic circuits, or filter the generated spin currents A^ 
One of the next most prominent goals of spintronics is to 
go further towards the microscopic regime^— and try 
to measure and manipulate single spins in quantum dots 
using spin biased circuits. Understanding the structure 
of spin current noise and response through quantum dots 
is therefore of primary importance. Moreover, the inter- 
play of strong interactions and the impact of quantum 
fluctuations of the spin on spin transport are also im- 
portant questions of fundamental interest on their own, 
and quantum dots, being the simplest strongly interact- 
ing systems, also play a prominent role in this regard: 
they allow for the systematic and controlled experimental 
and theoretical study of strongly interacting states. Al- 
though not easy to measure^"— dynamical noise spectra 
and response functions would allow to gain a clear insight 
to the structure of interactions.^^ 

In spite of its obvious importance, however, surpris- 
ingly little is known about the dynamical spin current 
response and the spin current noise spectrum of a quan- 
tum dot. Spin current shot noise in the sequential tun- 
nehng regime has been theoretically studied in Ref. |19j . 
and later in the co-tunneling regime by Kindermann^ 
However, these calculations focussed almost exclusively 
on d.c. properties, and have avoided the strong coupling 
(Kondo) regime^ which is much more difficult to reach 



and understand theoretically. 

In the present paper we make an attempt to character- 
ize the equilibrium spin noise spectrum and the dynami- 
cal spin response of a quantum dot. Here we focus on the 
local moment regime, ^^ where charge fluctuations can be 
neglected, and the dot can simply be described in terms 
of a local spin operator, S {S = 1/2), coupled to the left 
and right lead electrons through an exchange coupling, j. 
We shall study how time dependent spin polarized cur- 
rents can be injected through the dot at various temper- 
atures, and how these spin polarized currents fluctuate in 
time. To obtain a coherent and clear picture, we combine 
various numerical and analytical methods such as nu- 
merical renormalization group (NRG), perturbative and 
renormalization group calculations, and a master equa- 
tion approach. 

Though some of the results presented here are also 
valid in non-equilibrium, in this paper we focus exclu- 
sively on the equilibrium case, and leave the detailed pre- 
sentation of the rather technical non-equilibrium quan- 
tum Langevin calculation to a subsequent publication^ 
Nevertheless, even in this equilibrium case, the spin noise 
and the spin conductance display an extremely rich struc- 
ture; in addition to the temperature, T, and the Kondo 
temperature, T^, below which the dot spin gets screened, 
new time scales emerge. In the regime, T :^ Tk, e.g., we 
find that the Korringa rate (the characteristic spin decay 
rate). 
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plays a distinguished role. Furthermore, incorporating 
the effect of external spin relaxation processes also turns 
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FIG. 1: (Color online) Sketch of the setups to measure 
the spin-dependent left-right conductance, GYr (<^)5 and the 
noise, S^^ (uj). The arrows indicate the direction of the spin 
polarization in each lead. In the upper setup the spin filter 
on the right side detects the spin-resolved current at time t\ 
J^ (t^), induced by applying a spin-dependent voltage to the 
left lead at time t, VL^(t). The lower setup allows for inject- 
ing spins and measuring the spin-resolved currents in both 
electrodes. 



out to be important, and introduces an additional new 
rate, l/rg. 

In our analysis, we focus on two crucial quantities. On 
the one hand, we study the dynamical spin conductance 
G"" (co'). This quantity characterizes how a current J"(t) 
of carriers with spins polarized in direction n is generated 
in the left (r = L) or right (r = R) lead by a time 
dependent chemical potential shift, (5/i" (t) = eV^ (t), 
acting in lead r' on carriers polarized along n' (see Fig.[T]), 



{Jrit)) = I 



dt' C^!^, {t - 1') K" (f) ■ 



(2) 



Furthermore, we also investigate time dependent corre- 
lations of the currents J^(t), 

\t-t') ^ l{{j?{t),j:?'{t')}), (3) 
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and determine the corresponding noise spectra. Of 
course, in the equilibrium case studied here, G"" and 
6'"" are not independent, but are related by the fluctu- 
ation dissipation theorem [see Eq. (|26|) ]. 

One of our main results is, that - in the absence of ex- 
ternal spin relaxation and external magnetic field - both 
the noise and the conductance take on simple, universal 
forms, and apart from some geometry dependent prefac- 
tors, are characterized by just two universal functions. 
The left-right (r = L, r' = R) conductance, e.g., for 
n = az and n' = a^z reads 

Gli{uj) = j sm\^) [S^^rg{uj,T)^aa' g{uj,T)] , (4) 

where only the prefactor depends on the specific geome- 
try of the dot, characterized by the angle (j) (see Sec. HIl 
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— FIG. 2: (Color online) Sketch of the properties of the real part 



of the universal function, —g{uj,T) ^ GjJj^{uj), characterizing 
the cross-spin conductance through the dot. 



and Ref. 0), but the functions g and g are univer- 
sal functions of uj/Tk and T/Tk- The prefactor, -^ in 
Eq. (j4]), denotes the universal conductance quantum4^ 
The function g describes charge conductance through the 
dot, while g determines the cross-spin (a =t, a' =\) con- 
ductance, G]}^,{(jj). Our main goal is to study the proper- 
ties of g and those of the corresponding noise component 
in detail. The characteristic features of g are shown in 
Fig. [21 which gives a concise summary of our most im- 
portant results. As shown in Fig. [21 g vanishes in the 
limit, cj ^ 0, and \g\ develops a dip at cj < Tk-, and a 
broad resonance at a; ~ Tk for temperatures T <C Tk- 
The vanishing of the d.c. conductance is a consequence 
of the fact that spin can only be transferred from the 
spin-up channel to the spin-down channel by flipping the 
dot spin, S. Thus the amount of spin transfer is Hmited, 
and no d.c. spin-cross conductance is possible in the ab- 
sence of external spin relaxation. Temporarily, however, 
one can transfer spin between these two channels at the 
expense of accumulating spin on the dot. This amounts 
in the appearance of the broad resonance. 

The above-mentioned dip in \Re g\ also survives at 
temperatures, T ^ Tk'- There, by the simple argu- 
ment above, the cross-spin conductance also goes to zero 
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FIG. 3: (Color online) Real part of the conductance G[}j^{u) 
through the dot for T ^ Tk- The anomaly below the Ko- 
rringa relaxation rate, Ek, is a consequence of correlations 
between consecutive spin- flip processes, and is related to the 
dip in \^eg{uj)\ (see Fig.[2|). 
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FIG. 4: (Color online) Real time response of the current 
{JJi){t), as generated by a sudden change of amplitude, SVl. 



as cj -^ 0, however, this happens only below the Ko- 
rringa rate, oo < Ek [see Eq. ([!])], where correlations 
between consecutive spin- flip events become important. 
As a consequence of this, the spin-up - spin-up conduc- 
tance, G|^^(a;), develops a peak for frequencies uo < Ek- 
The relative size of the peak is universal (see Fig. [3]), and 
is determined by the ratio of spin-flip in the high tem- 
perature regime vs spin-diagonal scattering processes. 

By just analyzing the analytical structure of the con- 
ductance, G^^/ (cj), we can also make rather general con- 
clusions on the structure of transient response to a sud- 
den potential change at time t = 0, V^{t) = SV^ ©(^j, 
where Q{t) denotes the step function. The linear re- 
sponse is universal, and it depends only on tTx and 
T/Tk' Figure 2] sketches the structure of the zero- 
temperature response, {Jji{t)) generated by a sudden 
change of V^- The response is found to be logarithmi- 
cally singular at time t = 0, and at T = temperature, 
one observes at a time scale t ^ h/Tx an induced current 
bump of amplitude ~ {^'^ /h) SVl, followed by an expo- 
nential decay of the current response. The exponential 
decay we find is somewhat counter-intuitive: the long- 
time behavior is typically associated with Fermi liquid 
properties, ^^ and therefore one could naively expect an 
algebraic decay. However, the exponential decay found 
follows precisely from Fermi liquid theory, which implies 
certain analytical properties for the response functions. 

As mentioned before, the properties of the equilibrium 
noise spectrum are related to those of the conductance ji2£ 
Thus the noise has a structure similar to Eq. (j4]), and 
can also be described in terms of just two universal func- 
tions, s{uj) and s(cj), describing charge and cross-spin 
correlations, respectively. The left-right noise compo- 
nents S'£^ (cj), e.g., can be expressed as 



SrR{Lo) = -jTKsm^cP) [5, 



's{(jj^T) -\- aa s{(jO^T)] , 

(5) 

with the two universal functions s and s displaying 
markedly different behavior. At T -^ 0, e.g., s{uj) ~ 
\uj\/Tk, while spin cross-correlations scale as s{uj) ~ 
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FIG. 5: (Color online) Structure of the spin noise components 
5'J^(cj) and 6'l^(cj) for T^^Tk. 



This difference between s and s carries over to finite 
temperatures, where s remains finite in the d.c. limit, 
while the spin-dependent component s of the noise always 
scales to zero as cj ^ 0. In particular, at temperatures 
T ^ Tk the charge component of the noise is similar to 
what is found for an ordinary tunnel junction^*^ 

5(c.,T » Tk) « ^ j\u:,T) |- coth (^) , (6) 

and remains proportional to the temperature, T, for 
cj ^ 0. Here the only nontrivial physics is carried by 
the prefactor j2(a;,T) ^ l/ln^(max{|a;|, T}/Tk), which 
accounts for the Kondo-renormalized amplitude of the 
current, mediated by the exchange processes. Similar to 
^, however, the spin-dependent part s has a non-trivial 
structure, and exhibits a dip at frequencies below the Ko- 
rringa rate, uj < Ek^ where correlations between consec- 
utive spin-flip processes drive the cross-spin noise to zero. 
This amounts in somewhat unexpected low frequency 
anomalies in the noise spectra S^^, (a;), as sketched in 
Fig. [5] for the particular case of the left-right noise com- 
ponents. 

The results discussed so far hold in the case, where 
spin relaxation on the dot is induced exclusively by the 
exchange coupling, j, to the lead electrons. Then the 
d.c. cross-spin conductance and the shot noise vanish 
due to spin conservation. External spin relaxation, how- 
ever, slightly modifies the picture above, and in its pres- 
ence we find a small but finite spin shot noise and d.c. 
spin conductance. The precise dependence of the noise 
on the spin relaxation rate, I/ts, is rather complex, and 
is analyzed in Sec. IVII in detail. Experimentally, how- 
ever, the most relevant regime seems to be the one, where 
H/ts <^ T, Ek- There we find a residual cross-spin noise 
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For typical quantum dots, this residual noise and the 
corresponding conductance are thus very small compared 



to all other features, and therefore inclusion of l/r^ only 
slightly modifies the picture obtained within the Kondo 
model. 

The paper is organized as follows. In Sec. II we de- 
scribe the system's Hamiltonian and define the spin cur- 
rent operators. We then derive compact expressions for 
the spin noise and spin conductance in terms of universal 
functions using the Kubo formula. The zero temperature 
behavior is discussed in Sec. Ill, where we first show 
the analytical results in the perturbative and Fermi liq- 
uid regimes, and then present NRG results for universal 
functions obtained in the absence/presence of external 
magnetic field. As shown in Section III.B, application of 
an external magnetic field leads to a more complex be- 
havior and has somewhat similar effect as the external 
spin relaxation for spin currents polarized perpendicular 
to the applied field. 

Section IV is devoted to devoted to the study of the 
role of a finite temperature. In the perturbative regime, 
to determine the cj-dependence of universal functions we 
use a combination of master equation and perturbation 
theory, while in the strong coupling regime we resort to 
the Fermi liquid arguments. Transient response of the 
dot to a sudden switch of a time-dependent spin-resolved 
voltage is discussed in Sec. VI, while in Sec. VII we 
consider the effect of external spin relaxation on the dy- 
namics of spin correlations. Finally, the conclusions are 
given in Sec. VII. 



destroy electrons of spin a in the lead r. In terms of these 
fields, the Kondo interaction is given by 
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with j a dimensionless coupling related to the Kondo 
temperature, Tk ^ De~^^^, and a the Pauli spin ma- 
trices. Here the strength of the hybridization of the dot 
level with the leads is incorporated in the exchange cou- 
pling, j, while its left-right asymmetry appears through 
the dimensionless parameters, ^L/i?- These satisfy the 
relation v'j^ -\- vj^ = 1, and are parametrized by an an- 
gle (j) diS vl = cos(^/2), vr = sm{(j)/2).^^ Using this 
parametrization the maximal, T = temperature con- 
ductance is given by. Go = {2e^ /h) sin^ ((/)), where e^ /h 
is the conductance quantum. In the equilibrium case, 
studied here, the problem can be further simplified by 
noticing that only the combination ^ = vl "^l ^ vr iJjr 
occurs in Eq. (|9|). Therefore, rewriting the Hamiltonian 
in terms of the 'even' (^) and 'odd' (^) fields. 






COS ( - I iljLa + sin I - ) t/jR^ 



sin ( - 1 tl^La - COS ( - ) l/jRa 



(10) 



the interaction part simplifies to 



II. THEORETICAL FRAMEWORK 

A. Hamiltonian 

Throughout this paper, we focus our attention on the 
local moment regime of a quantum dot, where the dot is 
occupied by a single electron and it can be described in 
terms of a spin S, coupled to the leads. The Hamiltonian 
of the system consists of a part describing electrons in 
the leads, i^ieads, and an interacting part, i^int, which 
accounts for the exchange coupling between the dot and 
the leads, H = i^ieads + ^int- As usual, electrons in 
the leads are assumed to be non-interacting, and in the 
absence of a charge or spin bias, they are described by 
the Hamiltonian, 



Pleads = yZ ^ 4a (^) ^rai^) ds , 



(8) 



with cJcr(^) being the creation operator for a spin-a elec- 
tron with energy e in the left (r = L) or right (r = R) 
electrode and 2D being the bandwidth. The energy e 
is measured from the chemical potential, and the opera- 
tors Crai^) satisfy the usual anticommutation relations: 

The interaction part, Hmt, can be most easily con- 
structed in terms of the fields, ^ro- = /_£, Cra{s)ds, which 



ffi„t = I S (*to-*) 



(11) 



and the 'odd' field completely decouples from the dot. 



B. Spin current, spin conductance, and spin noise 

To calculate the spin noise and the spin conductance, 
we first need to define the spin current operators. Let 
us first focus on the case, where only the z-component is 
measured, and define the operators. 
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(12) 



measuring the total charge of the electrons in lead r and 
spin component, a, where e is the electron charge. The 
corresponding current operator can then be defined as 



j: 



dt 



(13) 



Notice that with this definition a positive current implies 
a charge flow towards the leads. Using the equation of 
motion we obtain 



j; 
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ra^ <^<^' H^r' a 
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h.c) . (14) 



The charge current is then simply expressed as the sum 
of spin currents, Jr = Xl^ J^ - 



At this point, it is useful to also express the current 
operators in terms of the even and odd fields, ^ and ^. 
Introducing the so-called composite fermion operators^ 
F(j = (Scr^)cr, we can write the current operator as a 
sum of two components. 
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(15) 



with the prefactors defined as '^l/r = (licos((/)))/4, and 
IL/R = ±sin((/))/4. 

It is instructive to analyze the properties of the current 
operators, Eq. ([15]). The components I^ satisfy 










(16) 



The second term in Eq. (p!9|) originates from the disconti- 
nuity of G^^/ (t) at t = 0, which can be shown to be real 
and just equal to (J^; J^, )u;=o- 

Let us now focus on the SU{2) symmetrical case. To 
simphfy the expression of G^^f , we can make use of the 
decomposition of /^, Eqs. ([Hj), and exploit the property 
(p!6|) as well as the fact that all cross-correlations of I^ 
and /^/ vanish. Taking furthermore the spin symmetry 
into account we find that G^^, reduces to the following 
form. 



G-;, (u;) 



h 



[5aa' arr' gi^J, T) + aa' ttrr' giui, T)] 



(21) 

Here the information on the geometry of the dot is ex- 
clusively carried by the matrices d{(j)) and a((/)). 



implying that the components I^ do not contribute to 
the charge current, Jr = "^Za^r- Furthermore, since the 
components /^, satisfy current conservation at the op- 
erator level, 



II + Ir 







(17) 



there is no charge accumulation on the dot, and Jhif) + 
Jnif) = is satisfied at the operator level at any time t. 
The current components, /^, on the other hand, do not 
satisfy current conservation, /£ + /^ 7^ 0, implying that 
spin accumulation is possible on the dot. These proper- 
ties follow naturally in the local moment (Kondo) limit, 
where charge fluctuations are completely suppressed, and 
only spin fluctuations are allowed on the dot. 
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while the functions g and g are completely independent 
of the dot geometry. They are defined in terms of the 
'reduced current operators'. 






z(^iF,-h.c.) 
z(^iF,-h.c.) 



and are given by the following expressions, 



/(^) = +^ {^TA^'h)^ - (2i-;2^)c^=o] , 



fl(w) = -^ (2^t;2^t)c^-(2t;2t)u,=o 



(22) 
(23) 



(24) 



C. Conductance and noise 

We shall now analyze the structure of the current re- 
sponse of the dot generated by an external perturbation, 
^R = V^{t) Qr{t)^ which shifts the potential of carriers 
with spin a in lead r. In linear response, the current is 
determined by the conductance G'^^, in Eq. (J2j), given by 
the Kubo formula, 



G;:,<(t,t') = -ie{t - t'){[j^{t),Qf,{t')]) , 



(18) 



This expression still contains the charge operators, which 
are non-local operators. To eliminate them, we differ- 
entiate with respect to time t' and use the equation of 
motion. After Fourier transformation we then obtain the 
following form of the Kubo formula. 



iuj G^^, {lo) — ( J^; J^, )oj — [Jr] Jr' )u;=0 , 



(19) 



where we introduced a compact notation for the Fourier 
transform of the retarded correlation function of any two 
operators. 
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(20) 



Here we have chosen a normalization such that for the or- 
dinary conductance through the dot we recover the well- 
known expression. 



GLRioo,T) = 2- sm\4,)g{u:,T), 



(25) 



with ^(cj, T ^ 0) = 1 in the unitary limit. 

Importantly, both g{uj^T) and giuo^T) are dimension- 
less functions, and, apart from some universal prefac- 
tors, they both determine physically measurable quanti- 
ties. By the basic principles of the renormalization group, 
this immediately implies that in the scaling limit, j ^ 0, 
D ^ 00^ and Tk finite, they must both reduce to func- 
tions of the ratios, uo/Tk and T/Tk- This can, of course, 
be checked explicitly by performing perturbation theory 
in the exchange coupling, j. 

In the equilibrium case, to which we restricted our- 
selves here, the symmetrized current noise [Eq. (J3j)] 
is simply related to the retarded Green's functions, 
{Jr]Jr')oji and thus to the conductance by the 
fluctuation-dissipation theorem. 



S^^> {uj) = -uj coth 



^)^eG^;;(c^). (26) 



Using Eq. (|2T]) . we then immediately find 



III. THE T = TEMPERATURE LIMIT 



where s and s are two real dimensionless universal func- 
tions characterizing the symmetrized equilibrium noise, 



s{uj,T) 
s{uj,T) 



The structures of the universal functions g, g, s, and s 
shall be thoroughly discussed in the following Sections. 



Before discussing the more complex case of finite tem- 
perature, let us focus on the much simpler T = case. 
There we can compute the universal scaling functions g 
and s "numerically exactly" using the machinery of nu- 
merical renormalization group (NRG),^^'^^ and can un- 
derstand all their important features relatively easily by 
combining perturbative renormalization group methods 
with Fermi liquid arguments. 



D. Non-collinearity 

The discussion above can readily be generalized to the 
case, where the polarizations of the injected and detected 
currents are arbitrary, but the system still exhibits SU(2) 
symmetry. In this case we first define the charge of car- 
riers polarized along n as 

nD -^ 

g^ = e V / cl^{e)-{l^na)^^^cra'{e)de. (28) 

/ J-D ^ 



The corresponding current operators, J^ = ^Q" can be 
expressed as 

J" = /" + /" ^ 

the even and odd current components being defined as 

I^ = ejjri{^^PnF-F^F^^), (29) 

with the projector Pn given by 

Pn = (l + nt7)/2. (30) 

Apart from the above changes of definition, the calcu- 
lations presented in Subsection III CI trivially generalize 
to this case, and the final expression of the conductance 
assumes the following simple form, 
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(31) 



with the functions g and g defined by Eqs. (|24|). Through 
the fiuctuation-dissipation theorem, we then obtain the 
following expression for the noise, 
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(32) 



with the functions s and s defined in Subsection III CI 
We emphasize again that the above expressions hold 
only in the presence of spin rotation invariance, and in 
an external magnetic field the conductance cannot be 
parametrized in terms of just two universal functions (see 
Sec. lmBl) . 



A. Analytical considerations 

Perturbation theory. The simplest one can do to com- 
pute the noise is to evaluate the current-current correla- 
tion functions order by order in j. The 0-th order con- 
tribution to Sl^ is found to be 
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g2 ^2j2 

— sin^ \uj\ -j^, (33) 



from which we can extract g and g using Eqs. ([2T]) and 
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(34) 



Evaluating higher order terms to S^^ , one obtains loga- 
rithmic corrections to g and g. The so-called leading log- 
arithmic corrections can be summed up by a perturbative 
renormalization group procedure^*^ which amounts in 
simply replacing j -^ j{uj) = l/ln(co'/Tx), and gives 
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, (35) 



for large frequencies, |a;| ::^ Tk- By the fiuct nation dis- 
sipation theorem, Eq. (f27|) . we then find for |a;| :^ Tk 
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(36) 



Fermi liquid regime. The above expressions approx- 
imate the scaling functions only at large frequencies, 
1^1 ^ Tk' In the opposite limit, |cc;| <C T^, perturbation 
theory in j breaks down. However, the relevant processes 
can be captured by a simple Fermi liquid approach^ In 
this very small frequency limit both g and g are analyt- 
ical. At T = temperature the charge conductance is 
unitary in the limit, cj ^ 0, and correspondingly. 



-LK 



(37) 



In contrast to ^, however, GIJ.,{uj -^ 0) and thus 
g{uj -^ 0) must vanish. This follows from the simple 
observation that it is impossible to generate a steady 
spin-down current in any of the leads by injecting only 



spin- up electrons. This simply follows from Fermi liquid 
theory^ According to this latter, at T = temperature 
the impurity spin is completely screened by the Kondo 
effect, and electrons injected at the Fermi energy are only 
subject to elastic scattering. As a consequence, a spin-up 
electron injected right at the Fermi energy conserves its 
spin. Electrons (quasiparticles), however, interact with 
each-other through a local interaction at the impurity 
site. For electrons injected with energy a;, this resid- 
ual electron-electron interaction generates spin-flip pro- 
cesses, and leads to a finite spin cross-conductance. By 
simple phase space arguments, the amplitude of these 
latter processes scales with the square of the energy of 



the incoming electron, ^ uj 
We thus conclude that 



and therefore g((jj) 



g = -a 



rn2 



s = —a 






(38) 



with a a universal parameter. 



B. NRG results 



In Section Hi B I we related the functions ^, and g to the 
correlation functions of the local operators, X^j and X^r, 
respectively [see Eq. ([M|) ] . At T = temperature, such 
local correlation functions can be accurately computed 
by NRG calculations^!^ which, as we show below, are 
indeed in full agreement with our previous analytical con- 
siderations. 



1. No external field, B — ^ 

Let us start by sketching how one can perform the 
NRG calculations in the case where we have no external 
magnetic field, B. To compute ^, we exploit the fact that 
correlation functions of the operators ^J-i^a and i^J^cr 
can be factorized onto correlation functions of F^ and 
^J.. The latter being trivial, after some algebra we then 
obtain for the real part of ^, 



3?e^(cj,T)= (39) 

? /^^^^^''^) [/(a;'-a;)-/(a;' + a;)] , 

with f{(jj) the Fermi function and qf{(jO^T) = 
— ^ ^m {F^;F1)^ the spectral function of the compos- 
ite fermion operator. The prefactor 7r^jf^/4 in Eq. (|39|) 
can be eliminated by observing that ^(0,0) = 1,^^ and 
therefore 



^'f 



0f{O,O) = 1 



(40) 



Using this relation, we can also express the real part of 
the scaling function g as 



^eg{uj,T) = 

2uj ^f(0,0) 



(41) 
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FIG. 6: (Color online) The zero temperature universal func- 
tions g^ g for the a. c. -conductance (upper panel) and for the 
spin noise s and s (lower panel) computed by NRG. 



with Qx^ x^i'^jT) = —^ ^rri {I^;ll)^ the spectral func- 
tion of the operator, Xa- = i(^J.Fcr — Fj^cr)- Computing 
thus the local spectral functions, qf and gx^x^ we can 
determine the real parts of the functions, g and g. From 
these, we can then directly determine the universal spin 
and charge current noise functions s and s. 

The real parts of g and g and the corresponding scal- 
ing functions, s and 5, as obtained by NRG^^ are shown 
in Fig. [6l The results presented were obtained by the 
Budapest Flexible-DMNRG code^^i^^ and clearly dis- 
play all features discussed in the previous subsection. All 
functions were plotted as a function of uo/Tk^ with Tk 
defined as a half- width of the function 3?e g. This lat- 
ter is essentially the a.c. charge conductance through 
the dot, studied in Ref. [36|, and displays the expected 
Kondo resonance, which appears as a peak at frequencies 

The behavior of the universal function, — 5Re ^ ~ ^lw 
is somewhat more surprising. It exhibits a maximum at 
a frequency, uo ~ T^, in agreement with the results of 
our analytical considerations, Eqs. ([35]) and ([38]) . This 
maximum corresponds to a temporary spin accumulation 
on the quantum dot at the "resonance" frequency, uj ^ 

The lower panel of Fig. [6] shows the universal functions 
for the spin noise on a logarithmic scale. The function s 
shows a clear linear behavior at frequencies uj < Tk- The 
deviations at larger frequencies from the linear behavior 
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FIG. 7: (Color online) Real parts of the universal functions gz and ^, and the corresponding noise scaling functions, Sz and s, 
as obtained by performing zero-temperature NRG calculations with a magnetic field Bz along the zth direction. 



are due to logarithmic corrections. At high frequencies 
—s follows the behavior of 5, and is linear apart from the 
aforementioned logarithmic corrections. For uo < Tk-, 
however, it deviates strongly, and scales to zero as ~ cj^. 



2. Effect of Zeeman field, B y^ 

The presence of a Zeeman field, Hb = —BS^ breaks 
the spin SU(2) symmetry. Therefore, for general spin 
polarizations, n and n\ the spin current noise and the 
spin conductance cannot be characterized by just two 
universal functions, as before. Using group theoretical 
arguments and exploiting the electron-hole symmetry of 
the Kondo model we can show that 5 universal functions 
need be used to characterize the complete n and n' de- 
pendence of C^^ . Moreover, these functions will not just 
be functions of uj/Tk and T/Tk^ but they also depend 
on the ratio, B /Tk- 

Here we do not venture to characterize all these func- 
tions. Rather, we focus our attention to the special case, 
where n \\ n' . Without loss of generality, we can then 
assume that n and n' are both parallel to the z axis, 
n = G z and n' = a' z. In this special case, one can 
show that, as a consequence of electron- hole symmetry 
and rotational invariance around the axis of the magnetic 



field, the contribution of the current I^ does not depend 
on the direction of the magnetic field, and that the total 
conductance has the structure, 

G^;, (c^, ^) = - {(5,,. ~arr' ~g{^, B) (42) 

+ cFcr' arr' [cos^(^) gz{uJ-, B) + sin^(^) g±{uj^ B)^ > , 

where is the angle between the magnetic field and n, 
and two new scaling functions replace the original scaling 
function, g{uj,T). Notice that all these scaling functions 
now depend on the size of the magnetic field, too. Appli- 
cation of the fluctuation dissipation theorem yields then 
a similar scaling form for the noise, with corresponding 
scaling functions, s, s^, and s±, 

S^^. (a;, 0) = -j TK{S^a' ~arr' 5(c^, B) (43) 

+ gg' arr' [cos^(6>) ^^(cj, B) + siv?{9) 5^(CJ, B)\ \ . 

To determine the scaling functions, g and gz^ and the 
corresponding noise scaling functions, we performed the 
zero-temperature NRG calculations for the case, where 
the magnetic field is parallel to the z axis, B^ = B (i.e., 
= 0). The results are summarized in Fig. [3 Applica- 
tion of a magnetic field gradually removes the spin degen- 
eracy of the dot, and leads to a splitting and suppression 
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FIG. 8: (Color online) Zero temperature NRG results for the 
real part of the universal function g± , and the corresponding 
noise scaling function s± for several values of magnetic field 
pointing along the x direction. 



a role similar to the external spin relaxation, discussed 
in Section IVIi and allows the impurity spin to flip back 
and forth, independently of the conduction electrons. As 
a consequence, the "pseudogap" feature in g±_ is gradu- 
ally filled up with increasing Bx-, reaching maximum for 



IV. THE FINITE TEMPERATURE CASE 

Although the NRG calculations presented in the previ- 
ous section could, in principle be extended to any finite 
temperature, a serious technical problem appears. As 
we mentioned in the introduction, the cross-spin conduc- 
tance, G£^^ and thus ^, must vanish at any finite tem- 
perature T in the cj ^ limit. However, currently used 
finite temperature NRG broadening schemes all produce 
linear (oc u) spectral weight in qx^x-^{oo^T)^ and there- 
fore by Eq. d^T]) lead to a finite and thus unphysical d.c. 
cross-spin conductance. 

We do not know of any way to get around this problem. 
Therefore, at finite T, we had to rely on analytical results, 
and combine them with the T = temperature results 
to obtain a coherent picture. At very high temperatures, 
T ^ Tk^ we can make use of perturbative approaches. 
However, as explained below, even in this regime simple- 
minded perturbation theory is insufficient, and we need 
to combine it with a master equation approach to ob- 
tain the complete uj dependence of the spin conductance 
and noise. Combining these perturbative results with 
scaling and Fermi liquid arguments, we are then able to 
understand the complete frequency and temperature de- 
pendence of the scahng functions, s, s, g and g. 



of the Kondo resonance for B ^ Tk- The splitting of the 
resonance can readily be observed in the scaling function 
g (i.e., the charge conductance through the dot). At very 
large fields, B ^ Tk, the conductance has a logarithmic 
tail for a; :^ B^ but is strongly reduced at frequencies 
uj < B^ where spin- flip processes are forbidden. 

Similar to g, the scaling function gz continues to vanish 
at cj = for any magnetic field. This property follows 
from the overall conservation of the total spin compo- 
nent parallel to the external magnetic field. The only 
effect of a magnetic field is thus to increase the region of 
suppressed conductance. 

Next, to determine the scaling function, g±^ we per- 
formed calculations with a field Bx applied along the x 
direction {6 = 7r/2). We also checked that g is the same 
as before, and does not depend on the direction of the 
applied field. The function g±^ however, is markedly dif- 
ferent from gz. Of course, for 5 = they are both equal 
to g. However, any finite magnetic field results in a fi- 
nite d.c. conductance, g±{uj = 0) 7^ 0. The reason of 
this is that in this case we measure a spin component 
which is not conserved for any finite magnetic field. In a 
sense, the perpendicular magnetic field in this case plays 



A. Perturbation theory and master equation 
approach 

At temperatures T :^ Tk corrections to the leading 
order perturbative results are small, and much of the 
noise spectrum can be understood based upon perturba- 
tive results. However, to understand the limitations of 
perturbation theory, we first need to understand the im- 
portant time scales in this high temperature limit and the 
way they influence spin transport. At T ^ Tk, trans- 
port through the dot occurs through individual exchange 
processes, whereby just one electron tunnels from one 
side of the dot to the other side of it. The typical time 
between such events is given by the "Korringa time", 
Tk, which we define as the inverse of the Korringa rate, 
Tk ~ H/Ek' To the leading order in j, it is given as, 
Tk ^ h/ pT. Tunneling events are, however, not instan- 
taneous in the sense that they are dressed by the in- 
ternal dynamics of the electron-hole excitations, created 
throughout the tunneling process. Correspondingly, the 
" duration" of a tunneling process is given by the thermal 
time, tt = h/T <^ tk- At very short times below the 
thermal time, t < h/T = tt (or at frequencies cj > T), 
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current-current correlations reflect just the internal and 
coherent dynamics of such a spin-flip event, as well cap- 
tured by the usual connected second order contribution 
to the current-current correlation function. 

As just stated, the usual bubble diagram only accounts 
for the structure of a single tunneling event. At times 
t > tk-, however, several independent incoherent tunnel- 
ing processes take place. These processes are correlated, 
since a spin-flip process that changes the dot spin from 
up into down (lt^4) must necessarily be followed by an 
opposite process when the dot spin flips from down to up 
(J|^^). These correlations turn out to be important for 
spin transport, and are obviously not captured by sim- 
ple perturbation theory. Fortunately, for times t > tt 
{(jj < T), the internal dynamics of a spin-flip event can 
be ignored, and we can make use of a master equation 
method as a complementary approach. There tunneling 
processes are taken to be instantaneous, just character- 
ized by some rates, but correlations between individual 
spin-flip events are properly accounted for through a clas- 
sical rate equation. 

From these simple arguments we thus conclude that 
the master equation approach must be valid for frequen- 
cies Co' < T, while simple-minded perturbation theory 
works for frequencies Ek < oj. Since Ek < T, the range 
of validity of these two approaches overlaps, as also con- 
flrmed by the actual calculations presented below and by 
the results of Ref. [H. 



1. Perturbation theory 

As explained before, for T :^ Tk and times t < tk 
(frequencies uo > Ek)-, perturbation theory (PT) is thus a 
good approximation, though it fails at longer times where 
already several spin-flip events occur, and the correla- 
tions between these spin-flip events cannot be neglected. 
Simplest 0-th order perturbation theory yields, e.g., for 
the left-right components of the symmetrized frequency- 
dependent spin noise 






e' . 2 . ^'j' 

sm Q uj coth 



uj coth ( — - ) + . . . , 
\2T) 



/ (jO > 

h ""' " 16 V2Ty 



with the dots referring to higher order corrections in j. 
The corresponding universal functions then read for uj > 
Ek, 






2 -2 



Tk 



coth 



.2t) 



n'f 



(44) 



Higher order terms give logarithmic corrections, and lead 
to a renormalization of j in these expressions. 

2. Master equation approach 

Let us now focus on the "classical" frequency regime, 
UJ <T. Here we can use a simple master equation (ME) 
approach, ^^ where we assume that, at any instance, the 
spin on the dot is either in a spin-up state 5 = ff or in a 
spin-down state 5* = ^, with corresponding probabilities, 
Vs{t). Conduction through the dot and spin relaxation 
are generated by scattering events, g, generated by the 
exchange interaction, Eq. (|9|). These scattering events 
are taken to be instantaneous, and consist of the scatter- 
ing of a spin a electron from lead r to a spin a' state in 
lead r' while changing the dot spin, S ^ S'^ 

q o {r\a\S' ^ r, a, 6'}. 



ir'a' -^r 



oc fv^vl,, fol- 



They occur with a rate, 7(g) 
lowing Fermi's golden rule. 

The dynamics of the dot spin is described by a simple 
master equation. 



dt In 



-r r 
r -r 



with the relaxation rate F 






^^4 



From Eq. 






given 
D) it 



as r 
follows 



(45) 



that 



spm polarization on the dot relaxes exponentially, 
{Sz) ~ e~^^^ Thus the rate F is related to the Korringa 
spin relaxation rate as, Ek.q = 2 F, which for the simple 
equilibrium case considered here takes on the usual 
expression. 



E 



K,0 



2F 



7T J 



(46) 



Here the label "0" indicates that this is just the leading 
order expression of the Korringa rate, and higher order 
terms in perturbation theory renormalize it [see Eq. ([T])]. 

To compute current-current correlations, we flrst no- 
tice that a scattering event q at some time r induces 
current pulses in the leads, J^{t) = AQ'^{q) S{t — r), 
with AQ^{q) G {±e, 0} the amount of charge transferred. 
Similarly, a series of events, {r^, qn} generates a current. 



J^it)=J2^QriQn)Sit-T„) 



(47) 



As a consequence, the stationary current-current corre- 
lation function can be simply expressed as 



{j;(t) Jf,m = 5{t) Y, P{q) AQUq)AQf,{q) + ^ P{q, t ; q' , 0) AQ-(q) AQ''J,{q') . 



q,q' 
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(48) 



Here the first terra describes the auto-correlation of indi- 
vidual scattering events, while the second term describes 
correlations between distinct tunneling events. The prob- 
ability P{q) in Eq. (|^8|) denotes the stationary rate for a 
given type of event, q = {r^cr^S' ^ r, <j, 6'}, and can 
be expressed as 



P{q) = l{q) Vs 



(49) 



with Vs the stationary probability distribution of the 
dot spin. The quantity P{q^t ; q'^0) denotes the joint 
probability rate of a scattering event q at time t and an 
event q^ at time t = 0. For the events q = {r2, cr2, 5*2 ^ 



ri,cri,5'i} and q' 
expressed as 



{r2,cr2,5'2 ^ ^i^'^i^'S'd it can be 



P{q,t; g^0) = 7(g) Ps.^S'S^) l{q') Vs[ , (50) 

where Vs^^s' {t) denotes the conditional probability that 
the dot spin evolves from state ^2 to Si during time 
t. The function Vs-^^s' (t) is determined by the mas- 
ter equation, Eq. (|4^ . and it is obviously this and only 
this quantity that generates time-dependent correlations 
between consecutive scattering events. Thus spin cur- 
rent correlations in this perturbative master equation ap- 
proach are directly related to the time evolution of the 
dot spin. 

Having set up this general framework, the detailed cal- 
culation of the classical noise spectrum is somewhat te- 
dious, but straightforward. Therefore, instead of pre- 
senting further details on it, let us just continue with the 
discussion of the final results. Within the master equa- 
tion approach, the left-right components of the spin noise 
read 
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From these equations we extract the fonowing approxi- 
mations for the universal scahng functions, 
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(52) 



(53) 



Remarkably, for T ::^ cj ::^ Ek,o, these results precisely 
coincide with the perturbative results, Eq. ([H]). This is 
indeed also clearly visible in Fig. [9l where we compare 
perturbation theory results for 6'£^ (uo) with results of 
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FIG. 9: (Color online) Upper panel: Frequency-dependent 
spin current cross-correlations {Sj^\^) and auto-correlations 
(^lr) i^ the absence of spin relaxation, as calculated using the 
master equation approach (ME, solid line) and second-order 
perturbation theory (PT, dashed line). In the calculations we 
assumed j = 0.2. Lower panel: Effect of a finite (rather large) 
external spin relaxation rate, I/ts = Ek,o/'2. 



the master equation approach. We remark that one can 
bridge these two approaches through a systematic but 
much more formal and difficult quantum Langevin ap- 
proach, already briefly sketched in Ref. |38|, a nd to be 
discussed in a subsequent publication, Ref. [23J . 

The fluctuation dissipation theorem, Eq. ([27|) . in this 
classical regime, uo <T^ amounts in the following expres- 
sions. 



^eg^^iu:) 



Tk 
2T 



, ME 
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UJ 



El 
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K,0 



3?e5^V) « ^S^^(c.) = ^j2. 



16 



(55) 



Notice that, in contrast to s and ^, the functions s and 
g are completely featureless in this frequency range. On 



12 



the other hand, in agreement with our ear her statement, 
Sj^j^{uj)^ 5Re Gj^j^{uj)^ s and 5Re g ah exhibit a dip be- 
low Ek, and scale to zero as cj -^ 0. There is a sim- 
ple heuristic picture behind this fact: a spin- flip process 
il^^ pumps spin- up electrons into the leads. However, 
it must necessarily be followed by a reverse spin-flip pro- 
cess, J|^^, where, on average, exactly the same amount 
of spin is pumped back as it has been pumped in before. 
These processes exactly balance each other in the long 
time limit, and lead to a vanishing cross-spin conduc- 
tance in equilibrium. We remark that if, however, there 
are external spin relaxation processes, then the spin may 
flip back spontaneously before pumping back the injected 
spin through the reverse spin-flip process. In this case, as 
we shall see in Sec. IVI) the conductance Gjjj^{uj) remains 
finite even in the uj ^ limit (see also Fig. [9]). 

A rather curious consequence of the dip in s is that, 
while S[jj^{uj) develops a dip below the Korringa rate, the 
noise component Sj^j^{uj) exhibits a peak of equal size, 
which precisely cancels the dip of —Sij^{uj) in the charge 
noise. This peak in SjJj^{uj) or the similar peak in Sj^j^{uj) 
may be more conveniently detected experimentally than 
cross-spin correlations. 



B. Beyond perturbation theory 

1. Logarithmic corrections 

In the previous subsection we discussed only the 
leading order perturbative and master equation re- 
sults. Performing, however, perturbation theory in 
j gives rise to logarithmic corrections. As long as 
max{T, IcjI} ;::^ T^, these corrections can be summed up 
using renormalization group methods ^i^^ and amount 
in the replacement of j by it's renormalized value, j -^ 
l/ln(max{T, \uj\}/Tk)- Apart from this substitution, 
however, the results of Subsections II V A ll and lTV A 2l con- 
tinue to be valid as long as T ::^ Tk- For uj > T^ e.g., 
we just recover the T = temperature results, Eqs. ([35]) 
and (j36j), while in the opposite limit, a; < T, we obtain 
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(56) 



with Ek = Ek{T) = 7iT/\i?{T/Tk) the renormalized 
Korringa rate of Eq. ([1]). Similarly, for the scaling func- 
tions g and g we obtain in this regime. 
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FIG. 10: (Color online) Sketch of the universal scaling func- 
tions s (continuous line), and s (dashed line), and the real 
parts of ^ (continuous line) and g{uj/TK^T /Tk) (dashed line) 
for T ^ Tk- The charge conductance and noise functions, 
g and s show no particular feature, while the spin scaling 
functions exhibit an anomaly below the Korringa rate, Ek- 



2. Fermi 



regime, T <^ Tk 



In the Fermi liquid regime, ^^ T <C T^, perturbation 
theory in j breaks down. However, we can derive the 
behavior of the scaling functions by two simple observa- 
tions. We first observe that in this Fermi liquid regime, 
j -^ oo^i^ and therefore the only remaining energy 
scales are T and Tk- Our second observation is that 
at the Fermi liquid fixed point, the residual electron- 
electron interactions are irrelevant^ and therefore phys- 
ical quantities are analytical functions of a;. In partic- 
ular, the asymptotic forms, Eqs. (|37|) and (|38|) remain 
valid up to the higher order corrections even at finite 
temperatures, 

?fie giu) = -a^ + 0{uj\ tW) , (58) 

^eg{Lo) = l + 0{uj^,T^) . (59) 

The scaling functions of the noise, s and s can then yet 
again be read out of the fluctuation dissipation theorem, 
Eq. dUl), yielding 



Fig. [TU] gives a concise summary of these results. 
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FIG. 11: (Color online) Sketch of the universal scaling func- 
tions s (continuous line) and s (dashed line) and the real parts 
of g (continuous line) and g (dashed line) for T <C Tk- 



These equations reduce to the T = expressions in the 
00 ^ T limit, and are valid as long as a; < Tk- Notice 
that in the cj ^ limit, the charge noise component 
s scales to a constant, oc T, while the spin component 
s scales quadratically to zero, just as in the T ^ Tk 
regime. Of course, for uo > Tk the scahng functions 
must also reduce to their T = temperature expres- 
sions, Eqs. (|35]) and (|36|) . The overall behavior of these 
functions for T <C Tk is sketched in Fig. [11] 
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FIG. 12: (Color online). Pole structure of the integrand in 
Eq. (|62p . for the spin-up - spin-down (upper panel) and spin- 
up - spin-up (lower panel) channels. In the spin-ti channel, 
the pole at cj = is canceled by the cj^ dependence of ^(cj, T), 
while in the case of spin-tt this pole survives and gives a finite 
response as t ^ oo. 



V. TRANSIENT RESPONSE 

Let us now turn to the discussion of real time tran- 
sient response, i.e., the time dependent current response 
in the left lead when a spin-dependent voltage of the 
form V]^{t) = 5V^ 6{t) is applied to the right electrode 
(see Fig. [1]). Within linear response theory, the average 
current pulse is just given by 
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with i/{uj — iS) the Fourier transform of the function. 
Surprisingly, just using Eq. ([62]) and the analytical 



properties of the functions G^j^{uj)^ we are able to 
make rather strong statements on the transient response, 
(^l(O)- Let us start by briefly summarizing these analyt- 
ical properties. First of all, being retarded response func- 
tions, GJ^(cc;), are analytical on the upper half plane. 
Moreover, as assured by Fermi liquid theory, they are 



also analytical in an extended region around cj = at 
any temperature. However, from perturbation theory we 
know that at very large frequencies, uj ^ Tk-,T^ they 
have logarithmic tails, \G^j^^{uj)\ ~ 1/\v?{uj/Tk) and 
thus tend to zero even in the universal scaling limit, 
I) ^ oo, Tk finite. Their asymptotic behavior and 
their symmetries [3?e G^^, {uj) = 3?e G^^, (~^) while 
^m G^^/ {uj) = — ^m G^^, {—uj)] imply the presence of a 
cut along the negative imaginary axis with an endpoint, 
-zA, with A (X max{T,TK} (see Fig. ^.^^ Further- 
more, as already discussed, G1^,{uj = 0) = 0, while the 
components G\l.,{uj = 0) = Grr'{^ = 0)/2 remain finite. 

Let us now discuss the properties of the response, 
Eq. (|62]) . First, we notice that due to the asymptotic 
l/\n{uj/TK) fall-off of Gl%{uj) and the analyticity on 
the upper half-plane, the integral contours in Eq. (|62]) 
can be closed upwards for any time t < 0. Therefore, 
(J£(t)) = for t < 0, i.e., it respects causality. The 
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response being zero even at t = is not entirely trivial: 
in the master equation approach, e.g., Gll,{uj) remains 
finite in the cu ^ oo limit, and one obtains an unphysi- 
cal jump at t = 0. We notice that the statement on the 
t = response being zero is equivalent to the Kramers- 
Kronig relation. Though the response (J£(t)) vanishes 
at time t = and is continuous for times t > 0, the slope 
of the response, ^(J£(t))|t=o is, however, infinite, since 
the integral J^ du Gll,{uj) logarithmically diverges. 

For times t > 0, the contours must be closed down- 
wards, as shown in Fig. [121 In the spin-up - spin-down 
channel, G]jj^{uj = 0) = 0, and therefore the pole at —i5 
does not give any contribution. The contribution of the 
cut to the spin up-down response can be written as 
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with Sg{x) = 2^m g{—ix + S) the cut of the universal 
conductance function, g{oj). Clearly, the contribution of 
the cut falls off as ~ e~^^ for long times. At T = tem- 
perature we have A ex T^, and furthermore Sg must be 
also a universal function, ^^(A -\- y) = Sg{y/TK)- There- 
fore, the response is a universal function of tTx- We 
can also tell the short time asymptotics of the response. 
Making use of the fact that the response is continuous at 
t = 0, we obtain for t <C 1/A the expression, 

p2 poo ^ 

{Jt{t))^-sm^<p 6g{A + y){e-y'-l)^. (64) 

1^ Jo y 

Since the cut scales for large energies as ^ 1/ln [y/TK)^ 
we get. 
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Remarkably, this result does not depend on the temper- 
ature, since it is determined only by the high frequency 
part of the conductance. Furthermore, since the length 
of the current pulse is determined by the exponential 
prefactor, r^ e~^*, we can also read out of Eq. ([65]) its 
height: for T <C Tk one has A ~ Tk-, and the height of 
the pulse is {Ji(t)) - W^. For T > T^, on the other 
hand, we have A ~ T, and the height of the current 
pulse, is ^ 5V]^/\v? ^. 

Figure [13] summarizes all the characteristic features of 
the current response (J^(t)), discussed above. The total 
charge pumped into spin-down channel of left lead is sim- 
ply given by the integral of the transient response, and 
is approximately 
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FIG. 13: (Color online). Transient current response in the 
spin-down (upper panel) and spin- up (lower panel) channels, 
upon a constant bias applied at t = in the spin-up channel. 



Remarkably, the coefficients appearing in this expression 
are just the high-temperature and low temperature ex- 
pressions of the spin susceptibility.^^ 

The analysis of the response of the spin- up carriers fol- 
lows very similar lines. The only major difference is that 
in this case the pole at —i5 gives a finite time independent 
contribution, and leads to an asymptotic response. 
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Otherwise, our discussions on the universal form of the 
response, and its short time 1/ln {t) singularity carry 
over to this case, too. Instead of giving further details 
on (j2(t ^ oo)), we just summarized its properties in 
Fig. [13 



VI. SPIN RELAXATION EFFECTS 

All results presented so far were obtained under the 
assumption that spin relaxation is generated by the ex- 
change coupling j, and there are no external sources of 
spin relaxation. In reality, however, external spin relax- 
ation channels are always present. In quantum dots, the 
dominant channel of (external) spin relaxation is usually 
due to hyperfine interaction between the confined elec- 
tron and nuclear spins in the host material, leading typ- 
ically to a dephasing time of the order of r^ ~ 10 ns or 
longer in the absence of magnetic field4Si^ These hyper- 
fine relaxation processes are thus characterized by an en- 
ergy scale hJTs '^1 — 10 mK, typically much smaller than 
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FIG. 14: (Color online) Shot noise Sl\^{uj = 0, T) as a func- 
tion of the spin relaxation rate -^ in the Fermi liquid regime 
(T ^ Tk) (upper panel) and in the perturbative regime 
lT:^Tk) (lower panel). 



the temperature. Coupling to piezoelectric phonons^ or 
electromagnetic fluctuations^ through spin-orbit inter- 
action or polaron dephasing processes due to coherent 
acoustic phonon generation are, in general, characterized 
by even longer dephasing times and smaller relaxation 
rates 4^ Therefore, for typical experimental parameters, 
we would naively expect l/r^ to be small compared to 
the temperature, T. Nevertheless, a finite r^ leads to 
qualitatively different results, since its presence lifts the 
constraint of spin conservation, and allows to have a fi- 
nite d.c. spin cross-conductance, G]jj^{uj = 0) 7^ 0. 

To investigate the effect of a finite r^, let us consider 
the perturbative (T ^ Tk) and Fermi liquid (T <^ Tk) 
regimes separately. In the regime T :^ Tk, we can read- 
ily extend our master equation analysis to include l/r^ 
and obtain. 
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Here we incorporated already logarithmic corrections by 
replacing the bare Korringa rate, Ek^q by its renormal- 
ized value, given by Eq. ([1]). Clearly, at large frequencies, 
uo ^ Ek, external spin relaxation does not play a role. 
However, for uj < Ek, it suppresses the dip in S[j^{uj), 
and leads to a finite cross-spin " shot noise" , 
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In other words, in the most relevant case, l/r^ < Ek, the 
cross-spin current noise is simply proportional to l/r^, 
as also found by Kindermann.^^ The full frequency spec- 
trum of the noise in the perturbative regime (T ^ Tk) 
is presented in Fig. [9l lower panel. 

In the Fermi liquid regime, T <C Tk, finite external 
spin relaxation also leads to a finite cross-spin noise. To 
estimate it, we first notice that I/ts just introduces a 
new frequency scale, and therefore we expect 



-^n 



G^LRi^ ^ 0) 



(69) 



K 



Notice that T does not appear in this equation. Corre- 
spondingly, the noise behaves for T < Tk as 
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The overall dependence of the cross-spin shot noise sig- 
nal on the spin relaxation rate is sketched in Fig. [T4j 
There we also display the physically not too relevant 
regime, l/r^ > T, where l/r^ becomes the dominant 
energy scale, and therefore we have j -^ l/ln(l/rsTx). 



VII. CONCLUSIONS 

In the present paper we studied the equilibrium spin 
current noise and spin conductance through a quantum 
dot in its Kondo (local moment) regime. We have shown 
that in the absence of external fields, they are both char- 
acterized by a pair of universal functions, and determined 
the properties of these functions. We have shown that - 
in contrast to the charge conductance [Glr] - the d.c. 
spin cross-conductance {G^^^) vanishes. Put in another 
way, there is no spin drag, and a spin-up current cannot 
generate a steady spin-down current, at least not within 
linear response. At T = temperature this obviously 
follows from Fermi liquid properties, and is thus valid for 
any interacting system with no spin-orbit coupling and 
with a Fermi liquid ground state. However, somewhat 
surprisingly, though it is not true for any interacting sys- 
tem, for a quantum dot, this property also carries over 
for finite temperatures. It is related to the simple struc- 
ture of the Kondo (or the underlying Anderson) mod- 
els, where spin transfer between spin- up and spin-down 
states can occur only through a single point, namely the 
dot state (or the dot spin in the Kondo model). There- 
fore, the spin currents generated by consecutive ^^^ and 
^^i\ flips of the dot spin precisely cancel each-other, and 
no d.c. cross-spin currents appear. Correspondingly, the 
cross-spin shot noise, SjJ-^{uj = 0,T) also vanishes at any 
temperature, and the noise spectra, S]^-^{uj) and S\J^{(jj) 
both exhibit related low frequency anomalies. 

External spin relaxation slightly changes the picture 
above. It partly removes the correlations between con- 
secutive spin- flip processes, and makes G^^(co' = 0) and 
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FIG. 15: (Color online) Correspondence between a single 
quantum dot device and a spinless double dot device. 



Sjjj^{uj = 0) finite. However, since the external spin-flip 
rate, l/r^ is typically much smaller than the other energy 
scales (T, T^, Ek), it only leads to small changes in the 
overall behavior of the noise and conductance functions. 

As we also demonstrated in detail, simple-minded per- 
turbation theory accounts only for the structure of in- 
dividual coherent processes, and fails badly to capture 
these correlations between consecutive processes, which 
happen to dominate the spin response at small frequen- 
cies. Therefore, one must be very careful when calcu- 
lating spin transport properties. Even in the pertur- 
bative regime, T ^ Tk^ simple-minded perturbation 
theory is valid only for frequencies above the Korringa 
rate, uo > Ek- To capture the physics at frequencies 
uj < Ek^ a supplementary master equation approach 
(valid for uo < T) can be employed. Alternatively, one 
can use a more systematic but also more technical quan- 
tum Langevin approach, which works for any frequency 
in the perturbative regime, T » T^, but neglects loga- 
rithmic corrections (see Ref. |23ll45| ). 

Although finite-frequency noise measurements are now 



available^— and spin polarized currents can also be rel- 
atively easily produced^^J^ measuring the low-frequency 
anomalies of spin cross-correlations, S'^}^(a;), seems to be 
a difficult task. However, the predicted low frequency 
anomalies are also present in the spin polarized conduc- 
tance, G|^^(a;), and noise, S\J^{(jj) (see Figs. [3] and [5]). 
These are experimentally much more easily accessible, 
since carriers must be polarized in the same direction. 
Alternatively, one can measure these cross-correlations 
in the charge sector^ by using capacitively coupled dou- 
ble dots (see Fig.[T5])4ii^ In the spin polarized case, the 
Hamiltonian of the double dot system maps to that of 
the Anderson model with anisotropic hybridization pa- 
rameters. Measuring noise or conductance between leads 
attached to the upper or lower leads of the double dot 
device shown in Fig. [15] is thus equivalent to cross-spin 
measurements in the single quantum dot setup. 
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